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Abstract 

We discuss activated escape from a metastable state of a system driven by 
a time-periodic force. We show that the escape probabilities can be changed 
very strongly even by a comparatively weak force. In a broad parameter range, 
the activation energy of escape depends linearly on the force amplitude. This 
dependence is described by the logarithmic susceptibility, which is analyzed 
theoretically and through analog and digital simulations. A closed-form ex- 
plicit expression for the escape rate of an overdamped Brownian particle is 
presented and shown to be in quantitative agreement with the simulations. 
We also describe experiments on a Brownian particle optically trapped in a 
double-well potential. A suitable periodic modulation of the optical intensity 
breaks the spatio-temporal symmetry of an otherwise spatially symmetric 
system. This has allowed us to localize a particle in one of the symmetric 
wells. 

Fluctuation-induced escape from a metastable state is at the root of many physi- 
cal phenomena, from diffusion in crystals to protein folding, and is closely related 
to nucleation in phase transitions and activated chemical reactions. In all these 
phenomena it would be advantageous to control the escape probability by ap- 
plying an external force. The problem of escape of driven systems has therefore 
attracted much attention in diverse contexts, a recent application being stochas- 
tic resonance We show that this problem can be solved in a very general 
form for a broad range of driving field frequencies, which goes far beyond the 
adiabatic limit. The analytic theory is compared with the results of analog and 
digital simulations. We then discuss experiments on controlling escape in mod- 
ulated optical traps. An important application of the results is the possibility 
of selective control of particles diffusion in a periodic potential, including both 
the rate and direction of the diffusion. 
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I. INTRODUCTION 



The question of how a system responds to an external field is one of the fundamental 
problems of physics. A strong nonlinear response is usually associated with a sharply reso- 
nant excitation of the system. However, the effect of external driving may also be extremely 
large for an important and wide class of phenomena related to large fluctuations, including 
escape from a metastable state and nucleation in phase transitions. 

The mechanism responsible is readily understood for adiabatically slow driving, where the 
driving frequency is small compared to the relaxation rate in the absence of fluctuations and 
the system remains in quasi-equilibrium. For systems in thermal equilibrium, the fluctuation 
probabilities are given by the activation law, W cx exp(— R/kgT). For large infrequent 
fluctuations, which are discussed in the present paper, the probabilities W are much less 
than all frequencies and relaxation rates. We will be specifically interested in activated 
escape, in which case R is the activation energy of escape. The driving force modulates the 
value of R quasi-statically and, even where the modulation amplitude \SR\ is small compared 
to R, it may still substantially exceed ksT, in which case W will be changed very strongly. 
We emphasize that the change of the activation energy is linear in the field amplitude, for 
\5R\ < R. 

For higher field frequencies, where the driving becomes nonadiabatic, the expected major 
effect of the field would be to "heat up" the system by changing its effective temperature. 
Indeed, in the weak- field limit, the escape rate W is known, theoretically and experi- 
mentally 0, to be incremented by a term proportional to the field intensity I rather than 
the amplitude A oc J 1 / 2 ||. However, one may ask what happens if the appropriately 
weighted field amplitude is not small compared to the fluctuation intensity (temperature), 
and whether an exponentially strong change of the escape rate will occur. 

Theoretical analysis of nonadiabatically driven systems is complicated, since one may no 
longer assume that the system is in thermal equilibrium. Whereas for equilibrium systems 
the exponent in the escape rate can be found, at least in principle, as the height of the 
free-energy barrier, for nonequilibrium systems there are no universal relations from which 
it can be obtained M; the situation with the prefactor is even more complicated [Q. Much 
effort has been put into solving the nonadiabatic response problem, in diverse contexts, and 
numerical results have been obtained for specific models (see e.g. |||7|). 

Recent theoretical results |||9[] show that, counter-intuitively, for high-frequency driving 
the change of R is proportional to the field amplitude, i.e., \nW is linear in A, over a broad 
range of A. The proportionality coefficient was called the logarithmic susceptibility (LS). 
Just like the conventional linear susceptibility, the LS relates the response of the system in 
the presence of external driving to its dynamics in thermal equilibrium in the absence of 
the driving field. We emphasize that the amplitude is a nonanalytic characteristic of the 
field, as it is obtained by taking the square root of the period-averaged squared field. We 
are therefore talking about a nonanalytic field dependence of the escape rate, and we need 
to determine a mechanism that would lead to such a dependence. 

In Sec. II we provide a general formulation which allows one to find, for a periodically 
driven system, the activation energy of escape induced by Gaussian noise with an arbitrary 
power spectrum. In Sec. Ill we outline the theory and analyze the frequency dispersion of the 
LS. We then discuss the results on the prefactor in the escape rate of a driven system |T(J and 
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analyze the full time- dependent as well as the time-averaged escape rate, including both the 
exponent and prefactor. In Sec. IV we present the results of analog and digital simulations 
of driven systems. These results provide a full qualitative and quantitative confirmation 
of the theory, and also reveal the underlying physics explicitly. In Sec. V we describe the 
experimental observations of the activated escape of particles in modulated optical traps. 
Sec. VI contains conclusions and a discussion of unsolved problems in activated escape, 
including the problem of statistical reconstruction of the dynamical model of a fluctuating 
system. 



II. GENERAL FORMULATION OF THE ESCAPE PROBLEM 

The idea underlying the theory of the LS is that, although the motion of the 
fluctuating system is random, in a large rare fluctuation from a metastable state to a remote 
state, or in a fluctuation resulting in escape, the system is most likely to move along a 



particular trajectory known as the optimal path (see |TT|-[19| and references therein). The 
effect of the driving field accumulates as the system moves along the corresponding optimal 
path, giving rise to a linear-in-the-field correction to the activation energy of escape. 

A natural theoretical approach to the escape problem is based on the path-integral 
technique. We will give a formulation which is based on this technique and allows one to 
find the logarithm of the escape rate for a periodically driven system. We consider a general 
case where fluctuations in the system are caused by a stationary colored Gaussian noise f(t) 
with a power spectrum $(u;) of arbitrary shape [18|,20 ]. The Langevin equation of motion 
is of the form: 

q = K(q;t) + f(t), K(q;t + r F ) = K{q;t), (1) 

where Tp is the period of the driving field. The noise is fully characterized by its correlation 
function <ft(t) = (f{t)f(0)) or by &(uj), the Fourier transform of <f)(t). The characteristic 
noise intensity is D = max$(o;)/2. 

If the noise is weak then, over the noise correlation time t corr and the characteristic 
relaxation time in the absence of noise t Te i, the system will approach the metastable periodic 
state q a {t) and will then perform small fluctuations about it |21[| . To escape from the basin 



of attraction of this state, the system should be subjected to a sufficiently large pulse of 
the force f{t). Various realizations of f(t) [the pulse shapes] can result in escape. Their 
probability densities are given by the functional p2[ 



V[f{t)\ = exp 



~ J dtdt'f(t)f(t-t')f(t') 



(2) 



where F{t) is a reciprocal of the noise correlation function (j)(t), J dt\ — — t') = 

D5(t - t'). For white noise, P{t) = <p(t)/4D = 5{t)/2. 

We assume that the noise intensity D contains a small constant, which is the small 
parameter of the theory. This parameter guarantees that the functional (0) is exponen- 
tially small for all pulses f(t) which can give rise to escape. In addition, its values differ 
exponentially for different appropriate f(t). Thus there exists a realization f(t) = / op t(£) 
which is exponentially more probable than the others. This optimal realization provides the 
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maximum to V subject to the constraint that the system ([I]) actually escapes. The path 
Qoptit) along which the system moves when driven by the optimal force f pt(t) is the optimal 
fluctuational path, q op t(t). 

From (H), the paths q opt , /opt provide the minimum to the functional 

1 r roo A poo 

n[q(t)J(t)] = - // dtdt' f( t )F{t-t')f{t')+ / dtX(t)[q-K(q;t)-f(t)}. (3) 

Z J J —oo J —oo 

They can be obtained from the corresponding variational equations of motion. The Lagrange 
multiplier X(t) relates f op t(t) and q op t{t) to each other. 

The boundary conditions for the escape problem follow from the fact that the system 
starts from the periodic attractor q a (t) = q a (t + t f ) in the distant past (on the order of 
the reciprocal escape rate), with / = asymptotically, and that, as the force decays after 
having driven the system away from the attractor, the system should not be brought back 
to the initially occupied basin of attraction. The latter condition is only satisfied if, 
for t — > oo, the system is approaching the unstable periodic state = q b (t + r F ) on the 
boundary of the basin of attraction to q a (t), 

f(t) -> 0, X(t) ^ for t -f ±oo; 

q(t) -> q a (t) for t -> -oo; q(t) -> q b (t) for t -> oo. (4) 

The time-averaged escape rate has the form 

W = C exp [-R/D] , R = min K. (5) 

The exponent R can be obtained for an arbitrary noise spectrum and an arbitrary periodic 
driving by solving the variational problem (||), (f|) numerically. In particular, in the case 
of white noise, where F(t) = S(t)/2, the Lagrange multiplier and the force f(t) can be 
easily eliminated from the variational equations, 2X(t) — f(t) = q — K, and the variational 
functional 1Z for the escape problem takes the form (cf. |§) 

1 r°° 

n[q(t)] = - dt[q-K{q-t)f. (6) 

4 J —oo 

The variational equations of motion for the problem (||]) are usually nonintegrable. In the 
case of a white-noise driven system this was pointed out by Graham and Tel H] . Generically 
there are several solutions which start from the attractor for t — ► — oo and arrive to a given 
state qj at a given time tf. The physically meaningful observable solution q op t(t) provides 
the absolute minimum to the functional 1Z |23| . 



The prefactor C in the escape rate (|5|) and the relation of W to a directly observable 
quantity, the time-periodic current from the basin of attraction, are discussed below. 

III. THE LOGARITHMIC SUSCEPTIBILITY 

We now turn to the case where the driving force F(t) is additive, 

K{q; t) = -U\q) + F(t), F(t + r F ) = F(t), (7) 
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and only weakly perturbs the system dynamics; in particular, it does not change the number 
of attractors or saddle states. Even in this case the effect of F(t) on the escape probability 



may be exponentially strong [pHTOl , because it is determined by the ratio of the field-induced 
increment 5R of the escape activation energy to small noise intensity D. We note that U(q) 
can be thought of as a metastable potential in which the system moves in the absence of 
periodic driving. 

To first order in F, the correction 5R can be obtained from the variational functional @ 
by evaluating the term oc F(t) along the zeroth-order path q^t(t), /opt(0> -^(t). However, 
special care has to be taken of the fact that the optimal escape path is an instanton ||. In 
particular, the function X^(t) is other than zero within a time interval of width ~ t C orr, ^rei 
and is exponentially small otherwise. At the same time, the optimal fluctuation leading to 
escape may occur at any time t c , in the absence of periodic driving (one can think of t c as 
the "center" of the instanton where |A^(t)| reaches its maximum). 

The field F(t) lifts the time degeneracy of escape paths. It synchronizes optimal escape 
trajectories, one per period, so as to minimize the activation energy of escape R. The 
field-induced change of R should be evaluated along such a trajectory, i.e. 

/OO 
dtx(t-t c )F(t) 
-OO 

= J2x(^F)F n exj>(inu F t c ), *(t) = -\ {0} {t), (8) 



where x(u>) = /f^, dtx(t) exp(iut), and F n is the nth Fourier component of the field [pOfl . A 
complete derivation for a white-noise driven system is discussed in Ref. |§; for the general 
case discussed here it will be given elsewhere. We note that Eq. (H) has a particularly simple 
form for sinusoidal driving, F(t) = Acosupt. In this case SR = — \x(u)p)\A. 

The change of the activation energy SR, and therefore the logarithm of the escape rate 
W, are linear in the field F(t). The coefficient x(u) is the logarithmic susceptibility (LS) 
||||. The function x is a characteristic of the system, as are, for example, the polarizability 
and other standard linear susceptibilities. It can be calculated for a given model or measured 
experimentally. 

Unlike the standard linear susceptibility which, by causality arguments, is given by a 
Fourier integral over time from to oo, x(w) is given by an integral from — oo to oo. The 
analytic properties of x{lv) therefore differ from those of the standard susceptibility and, 
in particular, their high-frequency asymptotics are qualitatively different. The standard 
susceptibility for a damped dynamical system decays as a power law for large w [e.g., as 
l/[U"(q a ) — iu], for the model (0)]. In contrast, from @ the LS decreases exponentially, 
x(uj) = Me~^ Tp , where r p = Imt p , and t p is the pole or the branching point of the function 
\(°>(t). The asymptotic behavior of x(u) is different, of course, if r p = 0, i.e. \(°'(t) has a 
singularity for real time. This happens, for example, if the potential U(q) has singularities 
encountered by the optimal path. Therefore it does not typically occur in dynamical systems. 

The LS takes a particularly simple form for a white- noise driven system. From @, (0), 

X (t) = -\®(t) = -q®{t), C = U'(C). (9) 

In this case, the explicit form of t p and the prefactor M in x{oj) are determined solely by 
the singularities of U'(q). They were obtained in Ref. 
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A. Complete nonadiabatic escape theory 



The notion of the LS makes it possible to find not only the exponent, but also the 
prefactor in the escape rate, and thus to obtain a complete nonadiabatic solution of the 
escape problem for dynamically weak driving. Since the celebrated Kramers paper |25| , the 



calculation of the prefactor has been one of the central problems in escape rate theory. For a 
periodically driven system, the escape rate W(t) is periodic in time. It can be introduced as 
a current j away from the metastable state, which is measured well behind the boundary q b 
of the attraction basin [for the model (^) with F = 0, q b is the position of the local maximum 
of the potential U(q)]. In the range |?7'(g)| F the current scales with q as 

j(q,t) = W[t-t d (q)\, dt d /dq = -l/U'(q). (10) 

Eq. (|10D provides a meaningful definition of both instantaneous and time-averaged escape 
rates. For weak driving, the values of the escape rate at different points q sufficiently far 
behind q b (t) differ only by a phase shift td{q), which makes it possible to make a sensible 
measurement of W(t) = W(t + 7>). For a white-noise driven system, an explicit expression 
for the time-dependent escape rate W(t) and for W was obtained [TO] by combining the 



results on the LS with the integral representation of the time-dependent probability density 
near q b . In particular, it was shown that 



r 2ir 

W/W = (2it)- 1 / d<f> exp[-6R(<f>/u F )/D], (11) 
Jo 

where Wq is Kramers' escape rate in the absence of modulation for an overdamped system 
(the type of systems which we discuss in this paper), and 5R(t c ) is given by Eq. (||). 

Since 5R(t c ) is a zero-mean periodic function, W always exceeds Wq. For small F/D, 
the correction to Wq is quadratic in F/D (cf. [0]). In the opposite limit of large F/D, the 
escape rate is changed exponentially, with ln[W'/Wo] ~ — D min SR(t c ), which coincides 
with Eqs. (|5|), (@). The dependence of the escape rate on time and the parameters of the 
system for a simple metastable potential is illustrated in Fig. [I]. 

The time dependence of the escape rate W(t) and the change of its form with varying 
parameters of the system, in particular with the frequency and amplitude of the driving 
force, were analyzed in Refs. |TD|,|2B[. The explicit form of the probability distribution in 
the vicinity of the boundary q b (t) was obtained in these papers, too. One of the conclusions 
which follows from the results is that the prefactor in the expression for the current j(t) 
calculated right on the boundary q b (t) has a totally different form from that in the current 
well behind q b (t), which gives the observable rate W(t) fllTf ). This is in contrast with what 
happens in the case of nondriven overdamped systems ||25|| . 

Calculating the current at the periodic boundary q b (t) was the goal of the recent papers 
by Lehmann et al. [27|. As noted before, the functional form of this current differs from 
that of the coordinate-independent instantaneous escape rate. In their analysis, Lehmann 
et al. adopted the idea [|5],|S| , Fig. ^|, of synchronization of optimal paths by a periodic field. 
The evaluation of the prefactor in |27| is based on an additional specific conjecture. Most of 
the specific results refer to a singular potential U(q) in Eq. (|7|): it consists of two opposite- 
sign parabolas matched between their extrema. However, the nonanalyticity of this potential 
should give rise to a deviation from the linear amplitude dependence of the activation energy 
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(0) for comparatively small amplitudes of the driving periodic force. The deviation will be 
strong where the amplitude of forced vibrations becomes comparable to the distance from 
the extrema of U(q) to the singular point where the parabolas are connected, as was indeed 
observed in Ref. |[27|| . However, as we showed earlier by solving the variational problem (^j), 
(|) exactly |5) (cf. Fig. 0), for generic analytic potentials the activation energy of escape is 
well described by the LS in a broad range of field amplitudes. We demonstrate this below 
by analog and digital simulations. 



IV. ANALOG AND DIGITAL SIMULATIONS 

A. Measuring the logarithmic susceptibility 

To test the relevance of the LS and to investigate its properties, we have built an analog 
electronic model [^| of the system (jl]) for the double-well Duffing potential 



U(q) = -\q 2 + (12) 

We drive it with zero-mean quasi-white Gaussian noise from a shift-register noise generator, 
digitize the response q(t), and analyze it with a digital data processor. We have also carried 
out a complementary digital simulation see Ref. |29f for details on the algorithm used and 



the noise generation. The analog and digital measurements involved noise intensities in the 
ranges D = 0.021 — 0.04 and D = 0.007 — 0.030 respectively, in dimensionless units. 

For escape from the state q a = — 1 of the white- noise driven Duffing oscillator, Eqs. (P), 



(|i2|) give the LS as 

X(oo) = TT^rKl - iw)/2]r[(2 + iu)/2)}, (13) 

where T(x) is the Gamma function. For sinusoidal driving, the measured time-averaged 
escape rate is compared with the expressions fly]), fll3|) in Fig. fj. We emphasize that the data 
refer to a strongly nonadiabatic driving, ujpt re \ = 0.6 (for the model (fP2"|), £ re i = 1/U"(q a ) = 
1/2), and cover the range from weak fields, A < D, to A/D = 10. The corresponding 
change of \8R\/D = \x{ojp)\A/D was < 4.2. The data and the theory are in full agreement, 
without any adjustable parameters. It is seen from the data that, for \x{ui F )\A/ D > 1 the 
dependence of In W on A becomes linear, as expected. We note that a qualitatively similar 
dependence of In W on the driving amplitude can be seen in the experimental data on driven 
Josephson junctions ||. 

In Fig. f| we show the data on the LS for several noise intensities. The activation energy 
R was obtained by measuring the slope of \nW vs 1/D. From @, the slope of R vs A 
yields the absolute value of the LS. The difference between the measured and calculated R 
arises from the noise intensity being not too small (D « 0.020 — 0.036 for the data points in 
Fig. |D), or in other words, comes from the field dependence of the prefactor in the expression 
for the escape rate W. As seen from Fig. ^|, when the latter is taken into account, there is 
full quantitative agreement between the theory and simulations. It is also seen from Fig. |3| 
that simulations with D as high as 0.03 still give the correct slope of SR vs A for large A, 
and thus the correct \x\. 
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The frequency dependence of |x(u;)|, a fundamental characteristic of the original equi- 
librium system, is compared with the theoretical prediction ( |13|) in the inset. As expected, 
the LS falls off exponentially at high frequencies, whereas the limit of x(oo) for uj —>■ cor- 
responds to adiabatic driving and can be obtained from the Kramers theory. We note that, 
generally, the LS is not a monotonic function of frequency: for underdamped systems, it 
displays resonant peaks ||. 



B. Switching between optimal paths 

We now turn to the investigation of a specific feature of the escape rate that is related 
to the minimization over t c in (H), It is expected to arise for a nonsinusoidal field, and in 
particular for a biharmonic one ||. Here, the periodic function SR(t c ) may have two minima 
per period. However, the activation energy will always correspond to the absolute minimum 
of 5R(t c ). For a certain relation between the parameters, the values of 5R(t c ) at the two 
minima are equal. The situation is then similar to the first-order phase transition where two 
minima of the free energy are equally deep. On the opposite sides of the phase transition 
line the system is in different states. In the present case, if the parameters pass through 
critical values where the minima of 5R(t c ) are equally deep, switching will occur from one 
minimum to the other. 

For biharmonic driving, a convenient control parameter is the phase difference 0i2 be- 
tween the field components Fi,F 2 . In the simulations we used the Duffing oscillator flT2|) 
driven by the field F(t) =0.1 cos(1.2t) + 0.3 cos(2.4t + (pu)- For such a field, the function 
6R{t c ) (|) has two minima. Their relative depths depend on (f>i2- 

The increment of the activation energy 5R = min 5R(t c ) as a function of 0i2 obtained 
from analog experiments and numerical simulations is compared to theoretical predictions 
in Fig. ^|a. For the critical value 0i 2 = cr , 5R has a cusp. On the opposite sides of the cusp 
it is determined by different minima of 5R(t c ). Relative numbers of escape events along the 
paths corresponding to these minima are shown in Fig. |5|b. The data clearly show that the 
contribution from one of the minima dominates everywhere except within a narrow vicinity 
of cr , where the contributions from the both minima are of the same order of magnitude. 

In Fig. |(c) we compare observed and predicted escape paths for 12 = cr (in the 
calculations, account was taken of the field- induced corrections). The coexistence of the two 
escape paths per period is clearly seen, and agreement with theory is excellent. 



V. DYNAMICAL SYMMETRY BREAKING IN A MODULATED BISTABLE 

OPTICAL TRAP 

A simple physical system which embodies fluctuation-induced escape is a mesoscopic 
particle suspended in a liquid and confined within a metastable potential well. The particle 
moves at random within the well until a large fluctuation propels it over an energy barrier. 
An optically transparent dielectric sphere can be readily trapped with a strongly focused 
laser beam, creating an optical gradient trap, i.e. "optical tweezers" p0| . Techniques based 



on optical tweezers have found broad applications in contactless manipulation of objects 
such as atoms, colloidal particles, and biological materials. Fluctuation-induced escape can 
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be studied using a dual optical trap generated by two closely spaced parallel light beams, 
as illustrated in Fig. ||. Such trap was implemented initially to study the synchronization of 
interwell transitions by low-frequency (adiabatic) sinusoidal forcing |nj . 



An important experiment with a particle in a double-well trap is a measurement of the 
transition rate in a stationary potential. Such an experiment can provide a rigorous test 
of the multidimensional Kramers rate theory with no adjustable parameters. Quantita- 
tive measurements require that the confining potential be adequately characterized. This 
can be done by measuring directly the full three-dimensional (3D) stationary probability 
distribution p(r) of a trapped Brownian particle [[3^ . 

A stable three-dimensional trap is produced by two focussed laser beams as a result of 
the electric field gradient forces exerted on a transparent dielectric spherical silica particle of 
diameter 2R = 0.6 /xm. Displaced typically by 0.25 to 0.45 /im, the beams create a double- 
well potential, with the stable positions of the particle centered at ri and r 2 . The stability 
perpendicular to the beam axis is due to the transverse beam profile gradient; in the beam 
direction the potential gradient is derived from the strong focusing of the objective lens [|30 



Relatively infrequent thermally activated random transitions between the potential wells 
occur through a saddle point at r s as depicted in Fig. ||. The experimental setup and the 
measurement technique have been discussed elsewhere fl32| . 

The full double-well confining potential U(r) is determined from the measured stationary 
distribution p(r) as U(r) = — k B T lnp(r). From the depths and curvatures of the potential 
wells and the curvature of U(r) at the saddle point r s , it is straightforward to calculate the 
Kramers escape rates. These rates can also be measured directly by placing the particle 
into one of the wells and measuring the average time it takes to switch to the other well. 
The potential U(r), and the barrier height in particular, can be systematically varied by 
changing the beam intensities. This results in an exponential change of the escape rate, 
thus making it possible to compare theory and experiment over a wide range of the escape 
rates. Extremely good agreement is obtained, as seen from Fig. [?]. 

The double-beam trap can also be used to investigate the effect of ac-modulation on 
transition rates. An interesting application of this effect is to direct the diffusion of a particle 
in a spatially periodic potential fl33| . It follows from the results of Sec. 2 that, for a generic 
periodic potential, the ac-induced change of the activation barrier differs depending on the 
direction (right or left, for example) in which the particle moves in escape. This makes the 
probabilities of transitions to the right and to the left exponentially different and results in 
diffusion in the direction of more frequent transitions. 

An effect closely related to directed diffusion, but more amenable to testing using optical 
trapping, is ac-field induced localization in one of the wells of a symmetric double-well 
potential. We expect both these effects to occur if the applied field breaks the spatio- 
temporal symmetry of the system [|J, |[34|| . The ratio of the stationary populations w±, w 2 of 
the wells is determined by the ratio of the period-averaged rates of the interwell i — > j 
transitions, 

wi/w 2 = W 21 /W 12 oc ex V ([5R 1 - 5R 2 ] /k B T), (14) 
where 5R\ 2 are field-induced corrections to the activation energies of escape from wells 1,2 



The experiment was conducted [35 for equal static barrier heights in the two wells AUi 
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AU2 = AU, with AUo set at ~ 7.5fcsT. The intensity of a laser beam was then modulated by 
an electro-optic device, giving rise to modulation of AU jksT with an amplitude ~ 2.5. The 
modulation frequency ujf/2tc was varied between 1 and 100 Hz, which covers the range from 
adiabatically slow to nonadiabatic modulation [the relaxation time is t Ie i ~ 10 -2 s]. This 
may be compared to the mean unmodulated transition rate Wq ~ 0.1 s _1 . Over this range, 
field-induced re-population was observed between the wells for a nonsinusoidal modulation 
waveform, so that ui\ ^ w 2 , 

The results on the instantaneous escape rates Wnif) and W2\{t) for an adiabatic mod- 
ulation (uf/Ztt = 1 Hz) are shown in Fig. || The barrier heights in the two wells were 
modulated in counter-phase. The form of the modulation was 5AU(t) = const x [sin(u;^t) + 
(1/2) smfiupt + 0i2)]- For this waveform, there is only one optimal escape path per period, 
for each well, and no switching between the paths occurred with varying 12 . 

As shown in the inset to Fig. [| the difference between the barrier heights in the two 
wells varies asymmetrically over the cycle. It depends on 12 and can be inverted if the 
phase angle is shifted by n. In other words, the modulated potential is not invariant under 
t — > t+ir/up, x — > —x, y — > —y (with x, y measured from the symmetry planes parallel to the 
beam axes, see Fig. |6|). It is this breaking of the spatio-temporal symmetry that leads to the 
escape rate from one of the wells being on average much bigger than from the other, as seen 
from Fig. || In turn, this leads to a higher population in one of the wells. Not only has the 
effect been observed for slow modulation, as evidenced by Fig. || but a population difference 
of 20% has been observed deeply in the nonadiabatic regime, with = 20 Hz, for the 

modulation amplitude used. This is sufficient to create significant directional diffusion, and 
demonstrates the onset of dynamical symmetry breaking. The dependence of W1/W2 on the 
phase shift 0i 2 is in agreement with the theory of Sec. Ill ||35|| . 

VI. CONCLUSIONS AND OPEN QUESTIONS 

We have shown theoretically, by analog and digital simulations, and by optical trapping 
experiments that fluctuations in driven systems, and in particular escape from a metastable 
state, can be effectively controlled by an external field. The field gives rise to a change 
of the activation energy of escape, which can be much bigger than the characteristic noise 
intensity (temperature) even for comparatively weak fields. Over a broad range of field 
amplitudes, this change is linear in the field even where the driving frequency exceeds the 
reciprocal relaxation time of the system and substantially exceeds the escape rate. The 
effect is described by a physically observable quantity, the logarithmic susceptibility. The 
LS relates the probability of large fluctuations in the presence of an external field to the 
dynamics in the absence of driving. It displays a specific frequency dispersion, which makes 
it possible to control, selectively, the escape rate of a targeted system. The LS can be 
calculated, for a given model of the system, and can be measured experimentally. 

An important application of the escape rate control is directed diffusion in a spatially 
periodic potential. As we have demonstrated, such control can be performed even for a 
symmetric potential, in which case not only the rate, but also the preferred direction of 
diffusion can be conveniently changed by changing the driving field parameters. 

The high efficiency of the escape rate modulation is largely due to the synchronization 
of optimal escape paths by the driving field. We have predicted and observed this synchro- 



10 



nization. We have also observed the related effect of switching between different branches 
of the activation energy as a function of the field parameters - a generic phase-transition 
type effect related to coexistence of different escape paths in systems away from thermal 
equilibrium. 

The results of this research are relevant to biological systems, since activated escape lies 
at the root of many biological processes at the molecular level, and modulation of the escape 
rate is often the way nature excercises control. However, detailed understanding of how this 
control is performed is missing in most cases. This is a fundamentally important and most 
challenging open scientific problem. 

Another important open problem of broad interest is whether large fluctuations can be 
used to learn about the dynamics of a fluctuating system away from stable or metastable 
states. The underlying idea here is that, in large fluctuations, the system explores remote 
areas of the space of its dynamical variables. An example where fluctuations have been 
used to find the global potential in which a system is moving was discussed in Sec. V. The 
problem becomes more complicated if the system dynamics is unknown. However, given 
that the system is most likely to move along a certain path during a large fluctuation, the 
observation of such paths can enable one to infer a dynamical model of the system. The 
process can then be iterated, as indicated by recent results |36| . 

The results described here are of interest also from the viewpoint of practical applications. 
An important example is the separation of colloidal particles and macromolecules. Our 
experiments show that selectively directed diffusion is a promising new approach to this 
problem. Another example is control of crystal growth using an ac field. The relevant 
nucleation rate will be changed by the field in a way similar to the escape rate. Therefore 
it should be possible to strongly modulate it both in time and space. 

The work was supported by the NSF through grants no. PHY-0071059 and DMR- 
9971537, and by the EPSRC (UK) under grants Nos. GR/L99562 and GR/R03631. 
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FIGURES 




A/D 



function of the scaled amplitude 
g 3 /3 [ 10 1 . The curves a to d refer to 



FIG. 1. The logarithm of the average escape rate ([Li] 
A/D of a sinusoidal field for the potential U{q) = q 2 /2 - 
the dimensionless frequency up = 0.1, 0.4, 0.7, 1.2. Inset: time dependence of the logarithm of the 
instantaneous escape rate for the same frequencies and A/D = 10 (0 = tupt), illustrating loss of 
synchronization of escape events with increasing ujp. 
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FIG. 2. Optimal escape paths (bold solid lines) of a periodically driven Brownian particle, 
i] = r\ — rf" + AcosLdt +f(t), for A = 0.1, w = 2 [from Ref. ||; 77 and w correspond to q and cj^t 
in the present paper, respectively]. The paths [given by Eqs. (||)— @] go from the stable to the 
unstable periodic states shown by bold dashed lines (by thin dashed lines, in the absence of driving). 
Thin solid lines show optimal paths in the absence of driving r](°\t — t c ) = — {l + exp[2(t — t c )]} -1 / 2 , 
with different t c . The driving lifts the degeneracy with respect to t c . The paths T)^(t — t c ) with 
the "right" t c [as given by (||)] are the ones around which the exact paths are oscillating. The 
linear nonadiabatic theory gives the decrement of the activation barrier to an accuracy 12%. 
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FIG. 3. The average escape rate W for a sinusoidally-driven Duffing oscillator (|T^) as a function 
of the field amplitude A, Wq is the escape rate for ^4 = 0. The driving frequency is ujf = 1-2, the 
white-noise intensity is D = 0.03. Solid line: the theoretical expression (|ll|); filled and empty circles 
are the data from analog and digital simulations, respectively, with no adjustable parameters. 
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FIG. 4. The dependence of the activation energy R on the amplitude A of the sinusoidal 
driving force with ojp = 1.2 for the Duffing oscillator as determined by electronic (open circles) 
and numerical (filled circles) simulations and Eq. ( |TT| ) (solid line) p4| . The data of analog and 
digital simulations refer to the noise intensities 0.028 < D < 0.036 and 0.020 < D < 0.028, 
respectively. The inset shows the absolute value of the LS of the system measured (open 

and filled squares for analog and numerical simulations, respectively) and calculated from (|l^) [full 
curve] as a function of frequency u. 
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FIG. 5. (a) The activation energy R as a function of phase difference (j)\2 with up = 1.2 for 
the Duffing oscillator driven by the biharmonic force F(t) = 0.1 cos(uiFt) + 0.3 cos(2uJFt + Mis- 
calculations based on (||) (full curve) are compared with data from electronic (open circles) and 
numerical (filled circles) simulations, (b) Relative numbers of escape trajectories following each 
escape path in the electronic (filled circles) and numerical (open circles) experiments compared to 
the calculated relative probabilities (full curve), (c) Measured escape trajectories for the electronic 
model (thin jagged lines) with the critical phase difference <pi2 = 4> C r ~ 3.57, compared to the 
calculated optimal paths (circles and triangles); solid lines are periodic states of ((!]), (0) in the 
absence of noise. The data was obtained with the noise intensity D = 0.028. 
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FIG. 6. Rendering of two focused laser beams, the equilibrium positions of the particle (rings), 
and a transitional path between the beams 




T T t meas / -1 N 

W (s ) 

FIG. 7. Comparison of the measured transition rates Vl7 meas and the rates calculated from the 
three-dimensional Kramers theory, W K , using the measured curvatures of the potential wells. The 
squares represent escapes from the well at ri and the triangles represent escapes from the well at 
Y2 in Fig. ||[ The line of slope one indicates the result expected if the three-dimensional Kramers 
theory correctly predicted the measured transition rates |32|| . 



18 




Phase of modulation (radians) 

FIG. 8. The least-squares fits to the experimentally determined instantaneous time-dependent 
switching probabilities for a particle in the adiabatically modulated double-beam trap, over a cycle 
ijj Ft of the modulating waveform. The phase angle between the first and second harmonics is 
012 = 7r/2. When the phase angle was incremented by ir, the escape rates from the left and right 
wells interchanged, within experimental error. Inset shows the instantaneous difference between 
the heights of the potential barriers in the two wells. 
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